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I .  INTRODUCTION 


An  important  objective  of  motion  sickness  research  is  to 
establish  quantitative  relationships  between  sickness  incidence  and 
various  parameters  of  the  physical  motions  that  induce  sickness.  An 
integral  part  of  this  objective  is  selecting  the  appropriate  motion 
parameters  upon  which  to  build  a  predictive  motion  sickness  incidence 
(MSI)  model.  Motions  conducive  to  motion  sickness,  however,  are  rather 
diverse  and  invariably  complex;  at  present,  only  whole-body  vertical 
sinusoidal  motion  has  been  studied  to  any  reasonable  degree. 

In  the  case  of  vertical  sinusoidal  motion,  research  has  linked 
MSI  to  the  frequency,  acceleration,  and  duration  of  exposure.  The 
apparent  relationship  between  these  variables  and  the  resultant  motion 
sickness  has  been  formally  described  by  an  empirically  based  predictive 
model  developed  by  Human  Factors  Research  (HFR) ,  Inc.  [3].  Using  the 
occurrence  of  frank  emesis  as  the  diagnostic  motion  sickness  criterion, 
the  HFR  model  predicts  the  cumulative  probability  of  emesis  as  a  function 
of  wave  frequency,  root-mean-square  (rms)  acceleration,  and  an  exposure 
time  of  not  more  than  two  hours  duration.  Initial  results  have  indicated 
that  the  effects  of  vertical  sinusoidal  motion  on  MSI  are  fairly  well 
described  by  this  model. 

Because  vehicular  motion,  be  it  on  land,  sea,  or  air,  tends  to  be 
broadband,  the  real-world  practicality  of  the  HFR  single  frequency 
model  is  limited.  Of  greater  practical  application  is  a  means  to  pre¬ 
dict  MSI  from  the  complex  motions  more  characteristic  of  operational 
environments.  Unfortunately,  the  current  paucity  of  data  regarding  the 


effects  of  complex  motion  spectra  on  motion  sickness  has  precluded 
the  development  of  a  predictive  MSI  model  for  general  motion;  it  is 
not  yet  known  what  characteristics  of  broadband  motion,  aside  from 
frequency  and  acceleration,  may  be  key  contributors  to  motion  sickness. 
Furthermore,  it  is  not  even  clear  for  broadband  motion  how  the  frequency 
and  acceleration  stimuli  might  best  be  quantified.  Further  research  is 
needed  to  resolve  this  perplexing,  but  important  problem. 

Although  it  is  customary  during  sinusoidal  motion  to  characterize 
the  dynamic  stimulus  in  terms  of  wave  frequency  and  rms  acceleration, 
other  equally  attractive  characterizations  do  exist.  In  this  paper  the 
predictive  utility  of  six  such  characterizations  is  examined.  Three  of 
these  representations  are  based  on  frequency  and  acceleration;  the  other 
three  are  based  on  frequency  and  jerk,  i.e.,  the  rate  of  change  of  accel¬ 
eration.  Although  each  representation  is  readily  generalizable  to  broad¬ 
band  motion,  none  completely  characterizes  the  more  complex  stimulus  pre¬ 
sent  in  broadband  motion.  Thus,  any  assessment  of  their  utility  for  pre¬ 
dicting  MSI  in  complex  motion  environments  will  be  tenuous  at  best. 

An  assessment  of  predictive  utility  is  further  complicated  by  the 
lack  of  complex  motion  data.  However,  to  obtain  some  preliminary  indica¬ 
tion  of  their  utility  for  predicting  MSI  in  broadband  motion,  an  experi¬ 
ment  involving  four  dual  frequency  waveforms  is  analyzed  in  this  report. 
Based  on  the  results  of  this  analysis,  possible  iiiplications  for  future 
research  and  model  construction  are  discussed. 


II.  SINUSOIDAL  DISPLACEMENT 


The  displacement  function  of  simple  harmonic  motion  is  given 
by  an  equation  such  as 

x(t)  “  Asin(ojt)  (2.1) 

where  A  denotes  the  wave  amplitude,  u>  »  2tt£ (uj  0)  denotes  the  angular 
frequency,  f  denotes  the  wave  frequency,  and  x(t)  denotes  the  displace¬ 
ment  at  time  t.  Equation  (2.1)  does  not  include  a  phase  constant  since 
such  a  constant  depends  only  on  the  initial  position  and  is  not  important 
for  the  description  of  sinusoidal  motion. 

The  acceleration  a(t)  corresponding  to  the  displacement  x(t)  is 
given  by 

a(t)  ■  x"(t)  *  -w2  As  in  (tot)  .  (2.2) 

In  sinusoidal  motion,  therefore,  the  acceleration  is  proportional  to  the 
displacement  and  in  the  opposite  direction.  Typical  parameters  for  measur¬ 
ing  acceleration  have  included  the  absolute  value  of  peak  acceleration 

(a  ),  the  root-mean-square  of  acceleration  in  each  half-wave  cycle  (a  ), 
max  rms 

and  the  time  integral  of  the  absolute  value  of  acceleration  imparted  in 

each  half-wave  cycle  (a)  .  Applying  these  definitions  to  (2.2)  yields 

a  ■  u>2  |  A I  ,  (2.3) 

max  11 

l 

3rms  “  I  1T/u)a2(t)dt]  2  -  /T5  uj2|a|  ,  (2.4) 

M  0 

and  . 

a  .  £  /  ja(t)|dt  -  (2/n)  u>2  |a|  .  (2.5) 

71  o 

It  is  clear  from  equations  (2.3),  (2.4),  and  (2.5)  that  a  predictive 
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rms 


MSI  model  based  on  a  ,  a  ,  or  a  could  just  as  easily  be  based , 

via  a  simple  transformation,  on  either  of  the  other  two  measures.  In 

this  sense,  then,  a  ,  a  ,  and  a  have  the  same  predictive  utility 
max  rms 

for  vertical  sinusoidal  motion. 

Consider  next  the  jerk  j(t)  (i.e.,  the  rate  of  change  of  acceler¬ 
ation)  corresponding  to  the  displacement  x(t) : 

j(t)  ■  a'(t)  -  -u)3Acos(tot)  .  (2.6) 

As  measures  of  jerk,  let  j  ,  J  ,  and  j  denote  the  analogs  of  a  , 
J  Jmax  rms  max 


a  ,  and  a,  respectively.  It  is  not  hard  to  show  that 


rms 


j  ■  u)3 1  A  |  , 

Jmax  1  1 

j  *  /73~ to3 1  A  I  , 
Jrms  1  1 


(2.7) 

(2.8) 


and 


j  -  (2/tr)  w3  |  A  |  .  (2.9) 

As  was  the  case  for  a  ,  a  ,  and  a,  it  is  evident  from  equations 

max  rms 

(2.7),  (2.8),  and  (2.9)  that  j  ,  j  ,  and  j  are  also  equivalent  in 

max  rms 

their  predictive  utility  for  vertical  sinusoidal  motion.  Furthermore,  it 

is  clear  that  conditional  on  frequency  the  predictive  utility  of  jmax  , 

j  ,  or  j  is  equivalent  to  that  of  a  ,  a  ,  or  a.  It  is  im- 
Jrms  n  max  rms 

material,  therefore,  which  of  these  six  motion  parameters  is  used  in 
conjunction  with  the  wave  frequency  f  as  a  basis  for  characterization 
or  prediction.  A  list  of  the  defining  relationships  in  which  all  the 
measures  are  equated  to  arms*3  presented  in  Table  1.  It  is  interesting 
to  note  that  for  sinusoidal  motion  jerk  can  more  simply  be  regarded  as 
the  interaction  between  frequency  and  acceleration. 

In  summary,  then,  the  researcher  has  at  his  disposal  a  number  of 
two-dimensional  characterizations  of  sinusoidal  motion  that  are  seemingly 
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Table  1.  Defining  Relationships  for 
Sinusoidal  Motion 
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different  but  essentially  equivalent  (up  to  a  change  of  units) ;  accord¬ 
ingly,  these  representations  offer  the  same  degree  of  predictive  utility 
for  vertical  sinusoidal  motion. 

The  HFR  single  frequency  model  can  be  used  as  a  predictive  basis 
in  actually  making  a  quantitative  assessment  of  the  common  predictive 
utility  of  a  ,  a  ,  a  ,  1  ,1  ,  and  j  .  For  example,  combi- 

nations  of  frequency  and  acceleration  for  which  data  is  available  can  be 
substituted  into  the  HFR  model  to  obtain  predicted  MSI  values.  These 
values  can  then  be  compared  to  the  corresponding  observed  MSI  values; 
the  apparent  degree  of  agreement  between  the  observed  and  predicted  MSI 
values  is  a  measure  of  predictive  utility.  Such  an  analysis  was  under¬ 
taken  in  [3]  where  a  6.1Z  rms  deviation  between  predicted  and  observed 
MSIs  for  a  number  of  data  conditions  was  observed.  This  represented, 
in  the  opinion  of  the  authors  of  that  paper,  "a  reasonably  good  fit  to 
the  data".  Their  analysis,  however,  was  made  on  the  same  data  that  went 
into  developing  the  HFR  model.  Hence,  the  apparent  degree  of  fit  is 
most  likely  an  exaggerated  indication  of  actual  fit.  A  truer  assessment 
of  predictive  utility  would  best  be  made  on  data  independent  of  that  used 
in  model  construction. 

After  discussing  vertical  sinusoidal  motion,  it  is  natural  to  in¬ 
vestigate  the  utility  that  the  various  parameters  may  have  for  predict¬ 
ing  MSI  in  motions  more  complex  than  sinusoidal  motion.  There  are  three, 
immediate  impediments  to  such  an  investigation: 

(1)  The  insufficient  data  base  for  complex  motion  precludes  a 
thorough  assessment  and  comparison  of  the  predictive  utility 
of  any  motion  parameters. 

(2)  Predictive  utility  is  contingent  on  other  motion  variables 
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present,  and  characterizing  the  pathogenic  effects  of 
broadband  motion  probably  requires  more  than  frequency 
and  acceleration  or  frequency  and  jerk  components. 

(3)  The  transformational  properties  that  played  a  key  role 

in  sinusoidal  motion  (cf.  Table  1)  do  not  hold  for  more 

complex  motion.  Thus,  a  ,  for  instance,  will  not  re- 

rms 

fleet  the  peak  (amav)  or  average  (a)  acceleration  during 

broadband  motion  in  precisely  the  same  way  that  it  does 
in  sinusoidal  motion. 

In  the  next  section  a  preliminary  analysis  of  some  existing  dual 
frequency  data  is  undertaken.  In  order  to  assess  predictive  utility, 
it  was  necessary  to  relate  dual  frequency  to  single  frequency  motion  and 
then  apply  the  HFR  single  wave  model.  Although  this  weakens  the  analysis, 
it  may  provide  some  preliminary  indication  of  predictive  utility.  At 
this  juncture,  there  seems  little  else  that  can  be  done  in  view  of  the 
above  impediments. 
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III.  DUAL  FREQUENCY  DISPLACEMENT 

The  displacement  function  y(t)  of  compound  harmonic  motion  is 
given  by  an  equation  such  as 

y(t)  -  A  sin(u)  t  +  6  )  +  A  sin  (u)  t  +  6  )  (3.1) 

1112  2  2 

where  A^  denotes  the  wave  amplitude,  -  2irf^  (  u)^  0)  denotes  the 
angular  frequency,  f ^  denotes  the  wave  frequency,  and  6^  denotes  the 
phase  constant  of  the  I*''1  (1-1,2)  component  sinusoid.  Without  loss  of 
generality,  it  can  be  assumed  that  one  of  the  phase  constants  can  be 
set  to  zero  and  therefore  omitted  from  equation  (3.1).  Hence  five 
parameters  are  needed  to  characterize  dual  frequency  displacement.  For 
predicting  motion  sickness,  however,  it  is  presently  unknown  what 
characteristics  (i.e.,  motion  parameters)  of  dual  frequency  motion  are 
important  predictors. 

Table  2  contains  data  from  a  pilot  study,  the  so-called  "Corre¬ 
lation  Study"  [1],  designed  to  determine  some  of  the  effects  on  human 
response  of  four  selected  dual  frequency  waveforms.  The  four  waveforms 
studied  were  created  by  summing  two  harmonically  related  sinusoids. 

The  resultant  waveforms  are  pictured  in  Figure  1. 

In  the  "Correlation  Study"  the  subjects  assigned  to  Condition  I 
(not  presented  here)  constituted  a  control  group  exposed  to  sinusoidal 
motion  at  the  fundamental  frequency  of  .17  Hz.  Condition  II  and  III 
added  a  second  harmonic  (.33  Hz),  differing  only  in  the  phasing  between 
the  fundamental  and  harmonic.  In  Condition  II,  0°  phasing  was  chosen  to 
minimize  the  peak  acceleration  of  the  resultant  waveform;  in  Condition 


MOTION  CONDITION 


(1) 


II 

III 

IV 

V 

' 

Fundamental  : 

Frequency  (Hz) 

.16 

.16 

.17 

.16 

rms  Acceleration  (g) 

.14 

.14 

.12 

.08 

Harmonic : 

Frequency  (Hz) 

.32 

.33 

.50 

.33 

rms  Acceleration  (g) 

.15 

.14 

.29 

.26 

Phase  Angle  (degrees) 

0 

+  90 

0 

0 

Number  of  Subjects: 

Total  Tested  ^ 

34 

34 

40 

40 

Quit 

2 

3 

9 

8 

Emesis 

16 

20 

21 

25 

No  Emesis 

16 

11 

10 

7 

MSI  (3) 

.53 

.68 

.75 

.83 

(1)  In  this  study.  Condition  I  was  a  single  frequency  control. 

(2)  The  duration  of  exposure  to  the  motion  was  120  minutes  unless 
vomiting  or  quitting  occurred  earlier. 

(3)  Quitters  included  as  part  of  the  emesis  group  in  calculating  MSI. 


Table  2. 


Results  of  "Correlation  Study":  Delivered  motion 
parameters;  number  of  subjects  who  yielded  data  or 
quit;  MSI  in  each  condition. 


ms:  Plots  of  displacement  versus 
(approximately  one  period) . 


Ill,  90°  phasing  was  chosen  Co  maximize  Che  peak  acceleration. 

CondiCion  IV  combined  Che  fundamenCal  frequency  of  .17  Hz  in  phase 
with  the  third  harmonic  (.50  Hz).  Condition  V  resembled  II  except 
that  the  level  of  acceleration  of  the  fundamental  was  reduced  while 
the  level  of  acceleration  of  the  harmonic  was  increased. 

Subjects  were  exposed  to  the  motion  for  120  minutes  unless  they 
quit  or  vomiced.  It  should  be  noted  that  a  supplemental  analysis  [5] 
to  the  "Correlation  Study"  has  studied  the  pattern  of  quits  and  con¬ 
cluded  that  "quitting  may  appropriately  be  regarded  as  part  of  the 
motion  sickness-prone  population".  Thus,  in  the  analysis  presented 
here,  quitters  were  included  as  part  of  the  emesis  group  in  calculating 
the  MSI  values  for  each  of  the  four  conditions.  In  effect,  it  is  assumed 
that  quitters  quit  because  emesis  was  imminent. 

In  order  to  relate  dual  frequency  to  single  frequency  motion 
it  is  necessary  to  reduce  dimensionality  from  five  to  two  dimensions. 

This  requires  trying  to  represent  the  spectral  content  of  dual  frequency 
motion  in  terms  of  a  frequency  and  either  an  acceleration  or  a  jerk  com¬ 
ponent.  To  represent  the  frequency  dimension  consider  the  mean  frequency, 
f,  as  determined  from  the  spectral  density  function  for  compound  harmonic 
motion: 

f  «  (A  2f  +A  2f  )/(A  2+A  2)  .  (3.2) 

112  2  1  2 

To  represent  acceleration  or  jerk,  a  ,  a  ,  a  ,  j  ,  j  ,  and  j 
can  be  calculated  directly  from  y(t).  These  six  parameters  together  with 
f  were  calculated  for  each  of  the  four  motion  conditions;  the  correspond¬ 
ing  values  are  presented  in  Table  3. 
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MOTION  CONDITION 


ii  m  IV  V 

Parameter: 

1  (Hz)  .170  .169  .186  .222 

a  (g)  .161  .151  .271  .238 

a  (g)  .201  .195  .314  .272 

rtns 

am  (g)  .354  .389  .506  .461 

max 

I  (g-Hz)  .278  .271  .844  .500 

1  (g-Hz)  .323  .321  .920  .545 

rras 

Jmax  (8"Hz)  '610  -  568  1Abu  *876 

Table  3.  Values  of  motion  parameters. 
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For  each  motion  condition,  the  values  of  f  and  each  of  the  three 
acceleration  and  jerk  parameters  were  substituted  into  the  HFR  single 
frequency  model,  suitably  reparameterized,  to  obtain  predicted  MSI 
values.  The  predicted  values  and  corresponding  p-values  are  given  in 
Table  4.  The  p-value  is  the  smallest  significance  level  for  which  one 
can  still  reject  the  hypothesis  that  a  given  parameter  combination  re¬ 
sults  in  the  prediction  of  a  true  MSI  value.  In  general,  the  smaller 
the  p-value  the  poorer  the  agreement  between  the  predicted  and  observed 
MSI  values. 

Inspection  of  the  p-values  reveals  that  of  the  six  parameter  com¬ 
binations,  Conditions  11  and  IV  are  fit  best  by  f  and  a  ,  Condition  III 
by  f  and  a  ,  and  Condition  V  by  f  and  j  .  The  latter  is  perhaps 
an  indication  of  a  shifting  mechanism.  Note  that  for  each  condition 
there  is  the  same  ordering  of  predictions  by  the  six  parameter  combi¬ 
nations. 

In  order  to  obtain  an  overall  indication  of  how  well  each  parameter 

combination  can  predict  MSI,  goodness-of-f it  tests  based  on  data  from 

all  four  conditions  were  performed.  The  resulting  p-values  were  as 

follows:  f  and  a  had  a  p-value  equal  to  .21,  f  and  a  had  a  p- 

value  equal  to  .15,  f  and  a  had  a  p-value  equal  to  .12,  f  and  Jmax  » 

f  and  j,  and  f  and  j  each  had  a  p-value  less  than  .001.  Overall, 

rms 

therefore,  f  and  a  had  the  most  success  at  predicting  MSI  in  the 
four  dual  frequency  motions.  Figure  2  presents  a  plot  of  the  p-values 
for  each  of  the  four  motion  conditions,  as  well  as  on  an  overall  basis. 


MOTION  CONDITION 


II 

III 

IV 

V 

Observed  MSI 

.53 

.68 

.75 

.83 

Predicted  MSI: 

f  and  a 

max 

.69 

(P-.038) 

.73 

(P-.515)* 

.80 

(p-,440) 

.74 

(P-.221) 

£  and  a 

.56 

(p-,692)* 

.54 

(P-.103) 

.75 

(P-.971)* 

.67 

(P-.034) 

f  and  a 

rms 

.61 

(P-.332) 

.60 

(P-.355) 

.76 

(P-.873) 

.68 

(p-,046) 

f  and  j 

max 

.84 

(P<.001) 

(P-.029) 

.96 

(p<.001> 

.83 

(P-.951) 

f  and  j 

.74 

(P-.005) 

.74 

(P-.430) 

.95 

(p<.001) 

.80 

(p«.694) 

f  and  j 

rms 

.76 

(P-.002) 

.76 

(P-.276) 

.94 

(p<.001) 

.79 

Cp-,637) 

*  Parameter  combination  producing  largest  p-value. 


Table  4.  Predicted  MSI  values. 


IV.  IMPLICATIONS  AND  FURTHER  DISCUSSION 


For  Conditions  II  and  IV,  the  acceleration  parameters  result  in 
reasonably  good  predictions,  while  the  jerk  parameters  do  not.  The 
opposite  is  true  for  Condition  V;  for  Condition  III,  both  types  of 
parameters  provide  roughly  equivalent  results.  It  would  appear,  then, 
that  both  acceleration  and  jerk  warrant  consideration  in  characterizing 
the  motion  sickness  mechanism.  It  should  be  stressed,  however,  that 
much  more  data  is  needed  in  order  to  explore  and  determine  fully  the 
relationship  between  dual  frequency  motion  and  the  resultant  motion  sick¬ 
ness. 

Unlike  other  approaches  which  have  been  proposed  [4],  the  use  of  f 

and  a  ,  discussed  herein,  results  in  creditable  predictions  that  can- 

rms 

not  be  rejected  at  conventional  significance  levels.  Nonetheless,  not 
enough  evidence  exists  to  conclude  with  a  high  degree  of  confidence  that 
the  HFR  single  frequency  model  can  indeed  be  satisfactorily  extended 
to  dual  frequency  motion. 

Although  the  approaches  examined  in  [4]  were  found  lacking,  it  should 
be  borne  in  mind  that  the  apparent  inability  of  a  weighting  scheme  to 
relate  dual  frequency  to  single  frequency  motion  for  the  purpose  of  ob¬ 
taining  MSI  predictions  could  imply  that: 

(a)  the  weighting  scheme  is  inappropriate, 

(b)  the  HFR  model  is  being  extrapolated  beyond  its  scope, 

or  (c)  the  relationship  between  dual  frequency  motion  and  the 
resultant  motion  sickness  is  simply  too  complex  to  be 
described  by  a  predictive  MSI  model  developed  solely 
from  single  frequency  data. 
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Because  of  data  limitations ,  the  results  of  this  paper  are  pre¬ 
liminary  in  nature.  With  the  aid  of  additional  data,  however,  a  more 
rigorous  investigation  of  predictive  utility  could  be  made.  Ideally,  a 
predictive  MSI  model  based  directly  on  broadband  motion  data  could  then 
be  developed.  Of  course.  Indirect  procedures  could  continue  to  be  refined 
and  examined  for  their  feasibility. 

In  either  case,  serious  consideration  must  be  given  to  the  design  of 
future  motion  sickness  experiments,  as  well  as  to  the  analysis  of  the  re¬ 
sulting  data.  From  an  analysis  and  model  development  standpoint,  it  may  be 
advantageous  to  place  more  emphasis  on  the  actual  time  to  emesis,  rather 
than  concentrating  only  on  the  occurence  or  nonoccurence  of  emesis  within 
a  given  period  of  time.  This  is  particularly  true  in  view  of  evidence  that 
time-to-emesis  may  be  regarded  as  a  random  variable  arising  from  a  mixture 
of  two  probability  distributions  [2J.  The  problem  of  predicting  MSI  can 
then  be  reduced  to  estimating  the  parameters  of  these  distributions. 

From  the  experimental  design  viewpoint,  not  only  must  the  motion  para¬ 
meters  of  major  interest  be  specified,  but  also  the  applicable  parameter 
space  must  be  well  covered.  Figure  3  illustrates  this  point  for  the  single 
frequency  case.  Consider  the  rectangular  region  illustrated  in  the  frequency 

and  a  space.  If  subjects  are  run  at  motion  conditions  corresponding  to 
rms 

the  four  vertices,  the  resulting  experiment  is  orthogonal  in  the  frequency 

and  a  parameters.  However,  if  the  region  is  mapped  into  the  frequency 
rms 

and  j  space,  the  rectangular  region  is  deformed  into  a  trapezoid,  re- 
rms 

suiting  in  the  loss  of  orthogonality  in  this  parameter  space.  This  might 

pose  serious  difficulties  in  model  construction  and  analysis,  particularly 

if  j  were  of  more  interest  than  a 
Jrms  rms 
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In  addition  to  the  problem  of  defining  the  most  appropriate  para¬ 
meters  (and  the  region  of  interest  within  the  parameter  space) ,  a  number 
of  trade-offs  must  be  made  in  the  quest  for  an  experimental  framework 
which  would  provide  the  most  useful  information  for  a  fixed  cost.  As  one 
example,  consider  two  experiments,  one  of  which  involves  many  experimental 
conditions  with  few  subjects  per  condition  and  the  other  which  involves 
only  a  few  experimental  conditions,  but  with  many  subjects  per  condition. 
In  this  situation  the  experimenter  is  faced  with  the  dilemma  of  balancing 
local  accuracy  against  global  coverage.  Unfortunately,  the  effects  of 
such  trade-offs  are  not  always  easily  quantifiable.  Research  on  the  im¬ 
pact  of  such  considerations  in  the  motion  sickness  prediction  framework 
is  currently  in  its  preliminary  stages. 
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